NEW  YORK  UNIVERSITY 

INSTITUTE  OF  MATHEMATICAL  SCIENCES 

LIBRARY 

4  WaAington  Mace,  New  York  3,  N.  Y 


INSTITUTE  OF  MATHEMATICAL  SCIENCES 

LIBRARY 

4  Waihington  Place,  New  York  3,  N.  Y 

NOTE  ON  THE  THEORY  OF  ROCKET  MOTORS 


As  a  basis  for  further  discussion  we  report  and  comment  on 
certain  results  in  a  paper  by  F.  Malina  concerning  the  theory 
of  rocket  motors*.   Schematically  a  rocket  motor  consists  of 
the  combustion  chamber,  in  which  a  gas  in  a  state  of  high 
pressure  is  created,  and  the  exhaust  nozzle,  through  which 
the  gas  flows  out  into  the  atmosphere.   Presently  we  are  only 
concerned  with  the  question  of  the  efficiency  of  the  exhaust 
nozzle.   The  main  problem  then  is  to  determine  the  propulsing 
thrust  resulting  from  the  gas  flow  for  nozzles  of  various 
shapes,  the  state  of  the  gas  in  the  combustion  chamber  being 
given. 

A  variety  of  simplifying  assumptions  are  made:  (i)  the  gas  is 
polytrophic  with  an  exponent y  =  1.2,  say;  (ii)  the  change  of 
state  of  the  streaming  gas  is  adiabatic,  losses  through  fric- 
tion at  the  walls  are  ignored;  (iii)  the  floor  through  the 
nozzle  is  steady;  (iv)  the  nozzle  diverges  so  gently  that  the 
classical  "hydraulic"  treatment  may  be  applied. 

A  main  result  on  the  basis  of  these  assumptions  is  that  the 
propulsive  thrust  is  a  maximum  when  the  floor  is  so  guided 
that  it  has  just  expanded  to  atmospheric  pressure  at  the  end 


*  Frank  J.  Malina,  "Characteristics  of  the  Rocket  Motor  Unit 
Based  on  the  Theory  of  Perfect  Gases." 
Journal  of  the  Franklin  Institute,  October,  1940 
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cross-section.   If  the  gas  is  expanded  below  atmospheric 
pressure  the  propulsive  thrust  decreases  very  rapidly. 

The  hydraulic  treatment  of  the  nozzle  flow  may  shortly  be 
described:  The  points  on  the  axis  of  the  nozzle  may  be 


characterized  by  their  distance  x  from  the  chamber  (c).   At 
each  cross-section  of  area  f  ,  the  flow  is  defined  by  average 
pressure  p  ,  density  y>   ,    and  velocity  v  .   The  flow  equations 
then  are: 

M  /Wx  =  m' 

the  mass  flux  per  unit  time,  is  independent  of  x;  so  is  the 


Bernoulli  constant 
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These  two  equations  together  with  the  adiabatic  law 


(3) 
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determine  three  of  these  quantities  as  functions  of  the  fourth, 
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When  the  pressure  p   is  taken  as  an  independnet  variable,  the 
graphs  of  v   and  f   are  as  in  the  diagram.   The  place  where 
the  cross-section  area  reaches  its  minimum  indicates  the 


"throat"  (t). 


F   =  (p   -  p.)  f   +  mv   ; 

X        X     *0     X        X  ' 


Let  F   denote  the  thrust  (in  positive  direction  opposite  to 

the  flow)  resulting  from  the  gas  pressure  inside  and  the 

atmospheric  pressure  p0  outside  exerted  against  the  chamber 

and  the  nozzle  up  to  the  cross-section  (x) .   Then  F   is 

x 

given  by 

U) 

i.e.  the  thrust  consists  of  a  combustion  due  to  pressure  dif- 
ference and  one  due  to  momentum  transport.   (Malina  refines 
this  formula  by  taking  into  account  that  the  flow  velocity 
will  be  radial  rather  than  axial  at  the  exit;  he  shows  then 
that  this  correction  is  so  small  even  up  to  nozzle  angles  up 

to  35   that  it  may  be  ignored.)   The  thrust  F   can  easily  be 

x 

expressed  as  a  function  of  the  pressure  p   or  the  cross-section 

x 

f  .   Malina  gives  the  result  in  dimensionl ess  form;  the  follow- 
ing formulation  deviates  slightly  from  his.   We  introduce  the 
quantities 
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x  = 
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then 


and 
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Malina  shows  the  graphs  of  the  thrust  ratio  Z  as  a  function  of 

the  cross-section  ratio  z  beyond  the  throat  for  various  valuees 

p0 
of  — .   Two  of  these  curves  are  indicated,  one  corres oonding  to 

Pc 


vanishing  outside  pressure,  the  other  to  a  small  pressure  ratio 

Pq 

—  =  .05.   One  observes  that  the  latter  graph  has  a  maximum  and 

Pc 

that  beyond  the  maximum  Z  drops  sharply  and  is  sensitive  to 

Po 
changes  in  — —  . 


The  maximum  of  Z  occurs  just  when  X  =  Xn,  i.e.  p   =  p  ,  as  fol- 


0- 


lows  e.g.  from 


dF   =  (p   -  p.)df   +  f  dp   +  mdv   =  (p   -  p  )df   . 
x    Vix    *Q'       x     x  *x       x     Fx    *0    x 

This  indicates  that  the  maximum  thrust  occurs  when  the  gas  at 
the  end  cross-section  has  just  expanded  to  atmospheric  pressure; 
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the  thrust  is  then  solely  due  to  momentum  transport.   (Malina 
does  not  mention  this  fact  explicitly.) 

Wnen  the  end  cross-section  is  so  large  that  the  exit  pressure 
would  be  below  outside  pressure  according  to  the  theory  as 
explained,  the  assumptions  on  which  this  theory  is  based  are 
no  longer  valid;  a  different  flow  pattern  involving  shocks  will 
occur.   Malina  determines  the  thrust  under  the  assumption  of  a 
straight  shock;  the  resulting  change  accounts  for  the  break  in 
the  indicated  graph. 

Malina  considers  as  an  example  a  case  with 

2c.  =  100,  T  =  11,100°F.   R  =  .0658  B.  T .  U  ./lb./°F. ,  y   =  1.2, 

and  finds  that  the  maximum  thrust  is  obtained  when  -g—  -   11.85, 

1t 

i.e.  when  the  exit  diameter  is  3.4-4  times  the  throat  diameter. 

The  maximum  thrust  is  given  by 

F 


Ve 


=  13.6, 


i.e.  the  thrust  is  13.6  times  the  force  which  the  atmosphere 
would  exert  against  the  end  cross-section.   Assuming  a  nozzle 
angle  0C   =  15°,  the  length  of  the  nozzle  would  be  9.1  times  the 
throat  radius  or  2.64-  times  the  radius  of  the  end  cross-section, 


It  seems  conceivable  that  operation  with  much  higher  pressure 

pc 

ratios  —  or  higher  temperatures  is  desired;  in  that  case  much 

p0 

longer  nozzles  would  be  required.   This  would  reduce  the  danger 
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of  over  expansion  and  shocks  within  the  nozzle   since  actual 
nozzles  then  would  rather  be  made  too  short.   The  hydraulic 
treatment  then  would  be  justified. 

The  design  of  a  nozzle  then  should  be  based  on  the  results  of 
this  treatment,  such  as  explained  in  Stodola's  book  on  turbines 
(cf.  section  30),  unless  other  physical  circumstances  are  to  be 
taken  into  account. 

When  a  nozzle  has  been  designed  in  the  indicated  way,  it  may  be 
desirable  to  determine  the  flow  through  this  nozzle  more  cor- 
rectly in  order  to  find  out  whether  or  not  the  hydraulic  treat- 
ment was  a  sufficient  approximation.   This  would  require  a 
mathematical  analysis  involving  the  solution  of  a  partial  dif- 
ferential equation,  which  would  be  rather  tedious  but  possible. 
It  seems  hopelessly  difficult,  however,  to  obtain  a  sufficiently 
large  manifold   of  such  solutions,  corresponding  to  various 
nozzles,  which  could  serve  as  a  basis  for  design. 


K.O.  Friedrichs 
October  13,  1943 


*  It  is  the  question  whether  or  not  the  jet  outside  of  the 
nozzle  is  to  be  considered;  in  that  case  the  ideas  of 
v.Neumann  about  the  flow  pattern  in  a  jet  should  be  taken  into 
account. 
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